IFSM-73-3A 



THE INCLUSION PROBLEM WITH 
A CRACK CROSSING THE BOUNDARY 


F. Erdogan 

AND 

G, D. Gupta 

ASE FILE 



Technical Report NASA TR-73-5 


March 1973 

National Aeronautics and Space Administration 

Grant NGR-39-007-011 





THE INCLUSION PROBLEM WITH A 
CRACK CROSSING THE BOUNDARY* 


F. Erdogan and G. D. Gupta 

Department of Mechanical Engineering and Mechanics 
Lehigh University, Bethlehem, Pennsylvania 18015 


Abstract 


The problem of an elastic plane containing an elastic 
inclusion is considered. It is assumed that both the plane and 
the inclusion contain a radial crack and the two cracks are 
collinear. The problem is formulated in terms of a system of 
singular integral equations. In the interesting limiting cases 
in which the crack tips approach the interface from either one 
or both sides, the dominant parts of the kernels become gener- 
alized Cauchy kernels giving rise to stress singularities of 
other than -1/2 power. For these unusual cases of a crack ter- 
minating at or crossing the interface stress intensity factors 
are defined and some detailed results are given for various 
crack-inclusion geometries and material combinations. 


1 . INTRODUCTION 

In studying the fracture of composite materials which con- 
sist of more than one perfectly bonded homogeneous elastic phase 
with different mechanical properties, it was shown that the 
singular behavior of the stress state in the close neighborhood 
of a crack tip does not remain "self-similar" as it enters and 
crosses an interface separating two phases of the composite 
[1,2]. If the crack tip remains in the same homogeneous medium 
as it propagates, during the crack propagation the characteristic 


*This work was supported by the National Aeronautics and Space 
Administration under the Grant NGR-39-007-01 1 and by The 
National Science Foundation under the Grant GK 11977. 


- 1 - 



square root singularity and the related angular distribution of 
the stresses at and around the crack tip remain unchanged, the 
only possible change taking place in the multiplicative constant 
known as the stress intensity factor. This makes it possible to 
apply any one of the conventional fracture criteria to this phase 
of the fracture propagation. On the other hand, since the singu- 
lar behavior of the stress field around the crack tip terminating 
at a bimaterial interface is drastically different than that of 
a crack tip imbedded into a homogeneous medium [1], as the crack 
enters and crosses the interface an abrupt change takes place in 
the crack tip stress field. Thus, since the stress field does 
not remain similar to itself during this phase of fracture propa- 
gation, for studying the related fracture phenomenon a closer 
examination of the crack tip stress field and some modifications 
of the existing theories or possibly a new fracture criterion 
are needed. 

A detailed treatment of this problem was given in [1] and 
[2] for the case in which the interface is a plane and the crack 
length and its distance to the interface are sufficiently small 
so that the perturbed stress field can be approximated by that 
of a crack in two bonded elastic half planes. However, in 
materials such as ceramics and fiber reinforced composites, the 
crack length is usually of the order of inclusion or fiber 
diameter and the stress state in the uncracked medium is quite 
different than that of two bonded half planes. Hence, for this 
type of problems clearly the assumption of bonded half planes 
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will not be valid. In this paper we will consider the plane 
elastostatic problem of a crack terminating at and crossing the 
bimaterial interface in an elastic matrix containing a circular 
elastic inclusion. The special case of the problem in which the 
crack is imbedded in the elastic matrix was recently discussed 
in [3]. Since [3] contains sufficiently detailed results of the 
single crack problem, in this paper, aside from a sample solution 
for the purpose of verification, we will not discuss this problem. 
Instead, we will give the solution of the problems of a crack in 
the inclusion with one or both ends approaching and terminating 
at the interface, of two collinear cracks one in the inclusion 
and one in the matrix, and of a crack crossing the interface. 

The analysis and the results given in [3] for the limiting case 
of the crack tip terminating at the inclusion boundary appear to 
be incorrect. Therefore some results for this case will also be 
given . 

2. THE INTEGRAL EQUATIONS FOR THE GENERAL PROBLEM 

Consider the plane elastostatic problem for an elastic 
matrix with constants k-j ,y-j containing a perfectly bonded cir- 
cular elastic inclusion of radius a and with constants 
where y i is the shear modulus, and = 3-4v n - for plane strain 
and k 1 - = (3-v.)/(l+v.j) for plane stress, being the Poisson’s 
ratio ( i = 1 , 2 ) . Let the medium contain two (radial) collinear 
cracks with end points at y = 0 and x = a^ , b-j , a 2 , t >2 (Figure 1). 

In addition to the geometry, let the external loads also be 
symmetric with respect to the plane of the cracks, y=0. The 
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figure 7. 


Inclusion-crack geometry. 


integral equations for this problem can easily be written down 
by using the dislocation solutions given in [4] and [5] as the 
Green's functions. In the usual manner the solution of the 
problem can be expressed as the sum of two sets of stresses: 

(a) the stresses in the medium without the cracks and under the 
given external loads, and (b) the perturbed stresses for the 
cracked medium where equal and opposite of the stresses found in 
(a) and applied to the crack surfaces are the only external loads. 
It is clear that the solution (b) alone will contain the singu- 
larities. 


Consider now the crack-inclusion problem shown in Figure 1. 
Let the crack surface tractions in the pertubation problem (b) be 


a lyy (x>0) = p l (x] 
a 2yy ( x ’°) = 


(a-| <x<b-j ) 
(a 2 <x<b 2 ) 


(1 .a,b) 


For example, for the uniaxial tension at infinity = a Q 

shown in the figure, the solution to problem (a), and hence, the 
tractions p-| and p 2 are given by 

a / „ \ r i a 


lyy 


(x,o) = a [l - 


a 2yy( X > 0 ) 


2x 
m ( k i + 1 ) 


2 »(K,-1) - (K 2 -l) a 4 3(m _ 1 ) 

,2 2m + ( k 2 - 1 ) ' 2 ( 1 + m < ] ) J » 


( 


2m + ic 2 - 1 1 + mi<i 


), 


p l ( X ) = - Q i yy ( x > ° ) » P 2 (x) = - °2yy ^ x * °^ » 

where m = y 2 / P i * Define 

f 1 ( X ) = gy [v 1 (x ,+0) - v^ (x ,-0) ] , (a 1 <x<b 1 ) 


( 2 . a - d ) 
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f 2 ( x ) = 97 [v 2 (x,+0) - v 2 (x,-0)] , (a 2 <x<b 2 ), (3.a,b) 

where and v 2 are the y-component of the displacement vectors 
in the matrix and in the inclusion, respectively, and f 2 can 
be looked upon as the unknown functions of the problem defined 
in the intervals (a^,b^) and (a 2 ,b 2 ), respectively. The consid- 
erations of displacement continuity require that 


k 

/ f k (x)dx = 0 , ( k=l , 2 ) . (4) 

a k 

Using the dislocation solutions given in [4] and [5] as the 
Green's functions, after some manipulations the following system 
of integral equations may easily be obtained to determine f-j 
and f 2 : 

b l f,(t) b l 

/ - 7.7 dt + / [ k n s ( x »t) + k llf (x,t)]f-,(t)dt 

a l a l 


b 2 tt(k,+1) 

+ y / [k 12s (x,t) + k ]2f (x,t)]f 2 (t)dt = — ^ Pl (x) , 

d £ 1 

( a -j < x < b -j ) , 

b 2 f (t) b l 

/ | - _7 ■ dt + 7 / [k 21 s (x,t) + k 21 f ( x , t ) ]f 1 (t)dt 

a 2 a l 


+ / [k 22 s (x,t) + k 22 f (x,t)]f 2 (t)dt = 


tt(ic 2 + 1 ) 

2 ~m> 


P 2 (x) 


(a 2 <x<b 2 ) (5.a,b) 


where 
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k ns^ x,t ) t-s ^ A 1 + A 2 ^ 2 x + 2 ( 3 s 2 - a 2 ) ( 1 - '^f) ] 

X 


t A ,[0-43) 1^4- ^V - -a l i 2 ] , 

1 ‘ x 2 (t - s) 2 a 4 t(t-s) 3 


k 1 lf(x>t) ^2 ^x + 2t tx 2 ) " ^ ( <2 + 1 ) ■ ] -* 2 ^ x 2 5 


A,+A„ -| „ A, -A. 2 2 

3 4\1 i /a a \ x , 3 4 X"d 


^ 12s ( x > t ) n - 2 } t-x + (A 3- A 4 ) trfcrr + 2 


t(t-x) c 


Ao - A* 

k- 12 f( x » t ) = —Jt + ^ A 5 


A ? + A 4 ] 3 2 
2 } x + ^ A 5 ' A 4^ 


A 4 ' A 3 a 2 ^ " A 3^ A 6 t 

2 x * 


2 tx 2 


A-l + B -1 y 1 

k 2 i s ( x > t ) C 1 " 2 + ^ A 1 " B l^t-*Fn 


V B 1 a 2 -x 2 


2 t(t-x) 2 


A 1 — B 1 -I 

: 21f (x,t) = [( 1 -m) M + g — L l^ 


Ao + A , 2 o c 2 9c -1 

k 22s^x» t ) = C 2 x ' a 3 2 ■ T^t^ 


+ *3 fl 0-T-> 


2 2 

4 sx 


A. 


s 3 (s 2 -a 2 ) 2 


( t- s ) 2 3 a 4 t(t-s) 3 


k fx t) = A3 + B3 + A r^- - a ^ t+3s > - 1 - A -M • 
K 22f u,t; 2x M 3 l 2 .3 t o 2 j • 

X L* X a 


( 6 . a-h ) 


y 2 1 +k-| 

Y ” vTj* l+< 2 * 


s = a /x , 


m = p 2 /y i , 


M = 


m( k-| + 1 ) 


A, = 


1 -m 


( i< 2 + m ) ( <2~ 1 + 2 m ) * 1 l+mK-j 


«2"niKi 

< 2+01 
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J 


1 


m-1 

m+K 2 5 

2 ( 1 - m ) 
Zm+K ^- 1 


mK, - k, 
A = 1 t 

m 4 1 +mK-| 


k^ + 1 


in( k-| +1 ) 


( 7 ) 


We note that for -a<a 2 , b 2 <a, a<a-j , (i.e., if none of the 
crack tips is on an interface) (see Figure 1) the kernels 
( k i j $ + k.j (i,j = l,2) are bounded functions of x and t in 
the intervals given by (5). Thus, in this case the set of equa- 
tions (5) is an ordinary system of singular integral equations 
with simple Cauchy type singular kernels. Since the displacement 
derivatives f-| and f 2 have integrable singularities at the end 
points of the corresponding intervals, the index of the equations 
is k=1 . Consequently, the general solution of the system will 
contain two arbitrary constants [6], which are determined from 
the continuity conditions (4). On the other hand, a close exam- 
ination of the "Fredholm kernels" k ^ ^ $ (i »j = 1 ,2) would 

indicate that, if one or more crack tips terminate at the bima- 
terial interface (i.e., if any one or a combination of the three 
cases a 2 = -a-j, b 2 = a, and a-j =a is valid) (see Figure 1), certain 
parts of these kernels become unbounded as both of the arguments 
x and t approach the end point on the interface. These parts of 
the kernels which go to infinity as (x,t)-»- + a are indicated by 
k---(x,t), (i,j = l,2) in (6). It is easy to see that k.. become 

I j S * J ^ 

infinite as (x-a) - ^, hence, together with the simple Cauchy 
kernels, (t-x)' 1 , they constitute a set of generalized Cauchy 
kernel s . 
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Singular behavior of the solution of integral equations 
having similar generalized Cauchy kernels was studied in detail 
in [1] and [2]. Following the complex function technique out- 
lined in [6] and using the procedure described in [1] and [2], 
if we define the unknown functions f-j and fp in terms of unknown 
bounded functions g-| and g 2 and unknown powers , 3-| > ' a 2 , and 
3 2 as follows: 

f 1 (x) = g-, (x) (b ] -x) 1 (x-a-j ) 1 , 

ou 

f 2 ( x ) - 92 ( x ) (^-x) ( x-ap ) » 

(-1 < Re(oij ,Bj ) < 0, j = 1 , 2 ) , ( 8 . a , b ) 

for various typical crack geometries, from (5) the characteristic 
equations giving and 6.. may be obtained as: 

(a) -a < a 2 < b 2 < a < a-j < b^ : 

cotira. = 0 , cot-rrB. = 0 , (j=l,2); (9) 

J J 

(b) -a < a 2 < bp < a , a = a-j < b-j : 

cotiTa-j = 0 , cotTrap = 0 , cotirBp = 0 > 

2costtB 1 + ( A-j + Ap) - 4A 1 (8 1 +1) 2 = 0 ; (lO.a-d) 

(c) -a < a 2 < bp = a < a-j < b^ : 

cot7m-| = 0 , c o t tt B - j = 0 , cot7rB 2 = 0 , 

2cosirap + (A^+A^) - AA^Cctp+l) 2 = 0 , (ll.a-d) 
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(d) -a = a 2 < = a < a-j < b-| : 

cotim-i = 0 , C0tTr3.| = 0 , = 3 2 5 

2cosTra 2 + (A 3 +A 4 ) - 4A 3 (a 2 + l) 2 = 0 , (12.a-d) 

(e) -a < a 2 < b 2 = a = a-j < b-j (through crack): 

cotuai = 0 , cotir3 2 = 0 , 3-| = a 2 , 

o TTOIq 9 2 

[(K-j - K 2 )cos — 2 — ” K-|(l+a 2 ) ] 

r\ TT Ot q pTTOtrt r\ 

+ K 3 sin cos' - K 2 (1 +a 2 ) = 0 , ( 1 3 . a - d ) 

1 + Kp 

K 1 = 2( 1 -m)/m , K 2 = — (1 + < ] ) , 

1 + KT 

K 3 - — JJ-2- + (1 + K, ) . (14) 

The equations (9), (lO.a-c), (ll.a-c), (12. a, b) and (13. a, b) 
give -0.5 as the acceptable root which is the well-known result 
for a crack tip surrounded by a homogeneous medium. (11. d) and 
(12. d) are identical to the previously obtained [1,2] character- 
istic equation for a crack tip terminating at a bimaterial inter- 
face. (13. d) is the same as the characteristic equation for two 
bonded quarter planes [2,7,8]. The characteristic equations 
(9 -13) are derived under the conditions that g.(a.j) and g ^ ( b - ) , 
(j=l ,2) are finite and nonzero. (13. d) is the expression of 
vanishing coefficient determinant in two homogeneous linear 
algebraic equations in g-j(a) and g 2 (a). Hence, in the case of 
a crack crossing the boundary g-|(a) and g 2 (a) are not independent 
and are related by [2] 
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g 1 (a) [Ai + A 2 - 4Ai (l+a 2 ) 2 + 2cosiTa 2 ] 
a-a 9 % 

= g 2 U)Y ( b~rf ) [a 3 + a 4 - 2 - 2 (i+a 2 ) (a 3 -a 4 )] . (15) 

The derivation of equations (9-15) follows very closely the 
procedure outlined in [1] and [2] in great detail and therefore 
is omitted in this paper. 


3. STRESS INTENSITY FACTORS 


From the viewpoint of applications of the results in frac- 
ture studies, one of the important quantities of interest is the 
strength of the stress singularity at the crack tips character- 
ized by the stress intensity factors. For the crack tips im- 
bedded in a homogeneous medium, the stress intensity factors are 
defined in terms of "cleavage" stresses and are related to the 
density functions f-| and f 2 as follows: 


k(bj) 

k ( a j > 


lim /2(x-b . ) a. (x,0) 

x-*b . J 

J 

lim /2 ( a . - x ) o. (x,0) 
x+a j J Jyy 


- lim /2 ( b . - x ) f • ( x ) to . , 

x->b • J J J 

J 

1 i m /2 ( x-a . ) f . ( x ) to . , 
x^a^. J J J 


w j = 2yj/(l+Kj), ( j = 1 , 2 ) . (16. a, b) 

The asymptotic expressions (16) may easily be obtained from (5) 
by noting that the expressions given by (5) are valid outside as 
well as along the cuts (a^,b-|) and ( a 2 > b 2 ) (i-e., 

p l(x) = a ]yy (x,0), ( -co< x <-a , a<x<°°) and p 2 (x) = a 2yy (x,0), 

(-a<x<a)), and by directly applying the f uncti on- theoreti c 
method to (5) [2]. Using the same procedure, from (5) the stress 
intensity factor for a crack tip terminating at the interface may 
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be obtained as follows [2]: 


-a < ag < bg < a = a-j < b-| 


k ( a ) = lim /2 (a-x) 
x->a 
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a 2yy (x ’ 0) 


“ ^2 9] (a) s i rttB-, 1+I-, ^ A 1 + A 2 - 2 - 2 ( 1 ) ( A ] -Ag ) ] > 


(17) 


•a < ag < bg = a < a-| < b-j : 


k(a) = lim /? (x-a) 
x+a 


-a. 


a lyy^ x *0) 


- ft g 2 U) sinTrag Kg+l ^ A 3 + A 4 " 2 " 2 ( 1 +a 2 ) ( A 3‘ A 4 > ^ 


(18) 


In the case of a through crack (i.e., if -a<ag<bg = a = a-j <b-| ) , 
for all practical combinations of material constants the func- 
tions f^ and fg, and consequently, the stresses have an integra- 
ble singularity at the intersection of the crack and the inter- 
face (i.e., -1 <ag = 8-| <0) . In fracture studies the quantities of 
interest here are the distribution of contact stresses along the 
interface. Thus, to characterize these stresses one may define 
the following stress intensity factors: 

V a > " yls y ° 2 °iyy (a - y) • k xy (a) = y 2 °ixy< a ->'> ■ 


(19. a, b) 

Developing asymptotic expressions for a. and a. around the 

j yy j xy 

point (a,0) and using (8), it can be shown that the constants 
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k and k are related to g-j(a) and g 2 (a) as follows [2]: 


V a) = 


P 1 r 9 l (a) m ' in . m -I 

^ 2 ^ 1 +mK, r "- L '-' ^ 


{■ 


7Ta o ' /c r 

2s i n— /b T a 


-I m+ic 2 


2 

g 2 (a) 

/a-a„ 




k xy< a > = 


"l 


g-| (a) 


Tra 


2cos 


2 

g 2 (a) 

/a-a„ 


[(l-ZcuWJ JL-] 

2 VE7T 2'1+mic, m+K 2 J 




(20. a ,b) 


4. EXAMPLES AND RESULTS 


Referring to Figure 1, if -a<a 2 <b 2 <a<a -j <b-| (including the 
special cases of single cracks, i.e., a 2 = b 2 or a -j = b ^ ) , the 
system of singular integral equations can be solved in a straight 
forward manner by using, for example, the technique described in 
[9]. In all the examples discussed in this section the external 
load was assumed to be the uniaxial tension O]” = o Q applied to 
the matrix perpendicular to the plane of the cracks and away from 
the inclusion-crack region (see equation 2). Following sample 
calculations for a simple crack were made as a spot check for the 
results given in [3]: 

(A) k-j = <2 = 2, P 2 /y-| = 1/3, a^/a = 1-1, b-j/a = 2.1 

k ( a 1 ) k(b 1 ) 

— = 1.482 , — = 1.1 60 , c -j = (b 1 -a 1 )/2 ; 
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(B) 


y 2 /J J i = V3, 


a , 


bj/a = 2 : 


(C) 


(D) 


k ^ - <2 2 j 

k (b, ) 

— = 1.233 , 


e 1 = -0.62090 
K-j = 1 • 8 , y 2 
k(b 7 ) 

— = 1.515 , 

k -j = 1.8, y 2 


k(a) 


-3- 


a o c l 


(see ( 1 0. d ) 


= 0, a -j / a 


1.092 , 

and (17)); 
= 1.05, 


c -j = ( b-j -a ) /2 , 


b -j / a = 2.05 : 


k ( a , ) 

— = 2.800 , 


c 1 ( b i “ a i )/2 ; 


0, a-| = a, b i =2 (edge crack) : 


k(b,) 

— = 2.808 , c, = ( b-, -a )/2 . 

1 1 

These results agree with that of [3] for a-j>a. However, because 
of the change in the power of the singularity 3-j for a-|=a, the 
extrapolated results in [3] are clearly in error (k(a-|) tends to 
zero or infinity as a-|-*a). 


For the cracks terminating at or going through the interface, 
the system of singular integral equations (5) (dominant parts of 
which have generalized Cauchy kernels) is solved by using the 
technique described in [10]. The results obtained for various 
material combinations and crack geometries are given in Tables 
1 - 6 and Figures 2-8 (see (16 - 19) for definitions of stress 
intensity factors). Table 1 shows the effect of on the 

power of stress singularity 3-j and on the stress intensity fac- 
tors for a crack in the matrix with one end touching the inter- 
face (the limiting case of the results given in [3]). Figure 2 
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1 1.5 x/a 2 

Figure 2. Crack surface displacement for a crack 
in the matrix with one tip on the interface 
(tc ] =k 2 = 1.8, b i / a = 2 , v = ( 1 +k ] )aa /u 1 ). 


Table 1. The effect of modulus ratio on the 

stress intensity factors for a crack 
terminating at the interface (a^=a, 

b i / a = 2 , k i =<2 = 1.8, Ci = (b-j-a)/2). 


U 2 

“ e i 

k(b 7 ) 

k ( a ) 



V,' 6 ’ 

0 


2.808 


0.05 

0.81730 

1.615 

1 .053 

1/3 

0.62049 

1.229 

0.5836 

1.0 

0.5 

1.000 

1.000 

3.0 

0.40074 

0.8610 

1.299 

10.0 

0.33277 

0.7969 

1 .389 

23.0 

0.30959 

0.7796 

1.375 

100 

0.29387 

0.7691 

1 .345 

300 

0.28883 

0.7667 

1 .348 
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shows the crack surface displacement v(x,0) for four typical 
values of ^2^^! ’ 0*05, 1 and 23) which is obtained from 

(see (3)) 

1 b] 

v(x,0) = - j / f-| { x ) dx , (21 ) 

where the normalizing factor is 

v 0 = ( 1+K i ) aa 0 /u l * (22) 


Tables 2-4 and Figures 3 and 4 show the results for a 
single crack located in the inclusion. The limiting stress 
intensity factors 0 and °° shown by an arrow in Table 2 is the 
trend based on the square root singularity. The correct stress 
intensity factors and the related or 3 2 are given in Table 3. 
Some of the results given in Table 2 are also shown in Figures 3 
and 4. The limiting values of the stress intensity factors 


shown in these figures for crack length 2 c 2 approaching zero are 
obtained from the uniformly loaded infinite plane solution with 
the stress state away from the crack given by (2.b), namely 


1 1m — - — 

c 2 - 0 a Q /^ 


1 _ 

2m+ic, 


(23) 


Table 4 shows the results for a completely cracked inclusion 


(i.e., a 2 =-a, b 2 =a). Table 5 and Figure 5 show the results for 
the case where both the matrix and the inclusion contain a crack 


(see the insert in Figure 5). The material constants used in 


this problem correspond to an epoxy matrix and an aluminum inclu- 
sion. In Table 5 the values of k(a-|) corresponding to a-|=a (the 
numbers in parentheses) are evaluated from (17) with g-j = -0.3381 1 
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0 0.5 c 2 /a 1.0 


Figure 3. Stress intensity factor for a symmetrically 
located crack in the inclusion ( K i = IC 2 = ^ ' ®) * 




Figure 4. Stress intensity factors for a crack located 
in the inclusion (k-| = < 2 = 1 . 8 , one tip fixed at 

= -0.9a or a„ = -0.75a, b 0 variable, c 9 = (b,-a,)/2) 


a„ 




k <t> 2 !/Vr 



factors for a crack in the 
in the inclusion (aluminum) 
= 23.077; a2=0.3a, b2 = 0. 


(b ] +a 1 )/2 


vari abl e) . 


Table 2. Stress intensity factors for a crack 

located in the inclusion ( <■■ = k 9 = 1 . 8 , 
c 2 = (b 2 -a 2 )/2 ). 1 


a 2 /a 

b 2 /a 

y 2 

M 1 

= 3 

y 2 

W 1 

= 1/3 

k ( a 2 ) 

k(b 2 ) 

k ( a 2 ) 

k(b 2 ) 


°o /S 2 


°o /F 2 

-0.9 

-0.75 

1 .324 

1.309 

0. 5886 

0.5950 

-0.9 

-0.5 

1 .451 

1.376 

0.5416 

0.5684 

-0.9 

-0.25 

1 . 572 

1.438 

0.5032 

0.5452 

-0.9 

0 

1.684 

1 . 501 

0.4719 

0.5219 

-0.9 

0.25 

1 .790 

1.572 

0.4450 

0.4969 

-0.9 

0.50 

1.890 

1.664 

0.4220 

0.4682 

-0.9 

0.75 

1.990 

1.822 

0.4020 

0.4300 

-0.9 

1 .00 

2.140 

00 

0.3830 

+ 0 

-0.75 

-0.5 

1 .314 

1.306 

0.591 7 

0.5950 

-0.75 

-0.25 

1 .389 

1.359 

0.5596 

0.5710 

-0.75 

0 

1 .475 

1.419 

0.5266 

0.5448 

-0.75 

0.25 

1 . 564 

1 .492 

0.4958 

0.5166 

1 

o 

CJ1 

0.5 

1 .655 

1 . 588 

0.4681 

0.4847 

-0.75 

0.75 

1.752 

1.752 

0.4437 

0.4437 

-0.75 

1.0 

1.907 

->■ 00 

0.4212 

-* 0 

-0.1 

0.1 

1 . 283 

1.283 

0.6046 

0.6046 

-0.25 

0.25 

1 .332 

1.332 

0.5796 

0.5796 

o 

LO 

o 

1 

0.50 

1.491 

1.491 

0.5144 

0.5144 

o 

VO 

o 

0.90 

2.062 

2.062 

0.3900 

0.3900 
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Table 3. Stress intensity factors for a crack 

located in the inclusion (k^ = < 2 = 1 . 8 , 
c 2 = ( b 2 -a 2 ) / ^ ) • 


a 2 /a 

b 2 /a 

-e 2 

~ a 2 

k ( a 2 ) 

k(b 2 ) 

-&2 

°o c 2 

-a 2 

a o c 2 



(Uj 

?/u-j ) = 3 



-0.75 

1.0 

0.5 

0.62049 

1 .907 

0.6175 

• 

o 

<o 

o 

1.0 

0.5 

0.62049 

2.140 

0.6300 

-1.0 

1.0 

0.62049 

0.62049 

0.7920 

0.7920 



(y 2 / 

'y-, ) = 1/3 



-0.75 

1 .0 

0.5 

0.40074 

0.4212 

0.9550 

• 

1 

o 

o 

1 .0 

0.5 

0.40074 

0.3830 

0.9400 

-1 .0 

1 .0 

0.40074 

0.40074 

0.9330 

0.9330 
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Table 4. Stress intensity factor for 

a completely cracked inclusion. 


y 2 

y l 

K 1 k 2 

= 1.8 

<1=2.2, 

<2 = 1 .8 

<1=1.8, 

<2=2.2 


= 2.2 

~ a 2 


-a 2 

k(a) 

" a 2 

k(a) 

” a 2 

k ( a) 

-02 

a 0 3 

-a 0 

°o a 

-a-> 
a o a 1 

- a 2 

a 0 a 

0.2 

0.36621 

0.7890 

0.38087 

0.7848 

0.32027 

1 .046 

0.33845 

1 .010 

0.6 

0.45025 

1.014 

0.47028 

0.9456 

0.42123 

1.174 

0.44466 

1 .068 

1 .0 

0.5 

1.0 

0.51991 

0.9209 

0.47724 

1 .107 

0.5 

1 .0 

no 

o 

0.57451 

0.8843 

0.59188 

0.8165 

0.55687 

0.9465 

0.57624 

0.8613 

5.0 

0.67885 

0.6555 

0.69124 

0.6194 

0.66380 

0.6940 

0.67733 

0.6500 
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Table 5. Stress intensity factors for a cracked inclusion 
in a cracked matrix ( y 2 / i = 23.077 , =1.6, 

<2 = 1.8, c-j - ( b-| -a-| ) / 2 , ^ 2 ~ ) • 


a l 

b i 

a 2 

b 2 

k < a i > 

k(b 1 ) 

k ( a 2 ) 

k ( b 2 ) 

a 

a 

a 

a 




a o /r z 

1.05 

1 .55 

0.45 

0.95 

0.335 

0.683 

1 .947 

2.716 

1.00 

1 .50 

0.45 

0.95 

(0.861 ) 

0.634 

1 .942 

2.732 

l'.OO 

2.00 

0.30 

0.80 

(1 .091) 

0.790 

1 .782 

2.061 

1 .25 

2.25 

0.30 

0.80 

0.681 

0.833 

1 .771 

1 .997 

1 . 50 

2.50 

0.30 

0.80 

0.831 

0.926 

1 .742 

1.919 

1 .75 

2.75 

0.30 

0.80 

0.898 

0.950 

1 .71 9 

1 .870 

2.00 

3.00 

0.30 

0.80 

0.932 

0.963 

1 .702 

1 .838 
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instead 


found from (10. d), (and with normalizing factor a Q c^ ^ 
o f C^/C^) • 

The results for the crack crossing the interface are shown 
in Table 6 and Figures 6-8. These results are also given for 
an epoxy matrix containing an aluminum inclusion. It should be 
noted that in solving the system of singular integral equations 
(5) for this problem, the single-valuedness conditions (4) are 
no longer valid. The two conditions necessary to account for 
the two arbitrary constants arising from the solution of the 
integral equations are the continuity condition v-|(a,0) = v 2 (a,0) 
and the relation (15) which must be satisfied by the functions 
g-| and g,,. The stress intensity factors k y and k xy given here 
are defined by (19) and are evaluated from (20). The limits 0 
and + °° shown by an arrow in the table (and indicated by dashed 
lines in the figures) toward which the stress intensity factors 
tend as the crack tip approaches the interface are again the 
consequence of the change in the power of singularity. For the 
materials under consideration the powers a. and 8-, ( j = 1 , 2 ) are 

J J 

found to be ( see (8) ) : 

-a<a 2 <b2=a = a-j<b^: a^=6 2 =-0.5, a 2 = 6-] =-0.27326 ; 

a l =a= * 3 2’ a 2 a < bi : a-j = -0.5, 8 2 ->-a 2 = 6i^-0. 33811 

3 1 ~ d = b 2 5 b-| a > a 2 : 3 2 = -0.5, a-j-»-6^=a 2 -* -0.33811 

(24. a-c) 

Figure 8 shows some sample results for the crack surface dis- 
placements obtained from (see (3)) 
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( k(Q),k(b))/V/c 


k(b ( )/V/c 


\ k(a)/(rc'^ 

xy ° 


b /a 


Figure 6. Stress intensity factors for a crack going 
through the matrix-inclusion interface ( k 1 - 1.6. 

k =1.8, Ug/M, = 23.077, «,= e 2 =-0.5, <» 2 = 6, = -0 . 27326 

c = (b,-a,)/2, a 2 = 0 fixed, b, variable). 






0.5 1.0 d / a 


Figure 7. Stress intensity factors for a crack 
going through the interface (^=1.6, k 2 =1 - 8, 

y 2 / Pl = 23.077, a 2 = 0 1 = -0.27326, 2c = (b 1 -a 2 )=a 

fixed, d= (b- | +a 2 )/2 variable). 


Table 6. The stress intensity factors for a crack 

crossing the interface (k^ =1.6, = 1 . 8 , 

y 2 /y 1 = 23.077, a 2 = -0. 27326, c = (b ] -a ] )/2 ). 


a 2 

■a 

k ( b i ) 

k ( a 2 ) 

1 

mu 

a 

a 

0 o /f 


<V'“ 2 

-a 2 
o Q c * 

0 

1.0 

1 


m 

00 

0 

1.1 


3.564 


0.170 

0 

1.3 


3.701 

BSE9 

0.0894 

0 

1.5 

0.570 

3.756 


0.0565 

0 

1.7 

0.626 

3.799 

- 0.171 

0.0342 

0 

1.9 

0.672 

3.838 

" - 0.0835 

0.0167 

0 

2.1 

0.710 

3.874 

- 0.0113 

0.00227 

0 

2.5 

0.765 

3.935 

0.105 

- 0.0211 

0 

3.0 

0.811 

3.996 

0.219 

- 0.0440 

- 1.0 

1.5 


oo 

-> - 00 

-> OO 

i 

o 

10 

1.5 

0.920 

5.600 

- 1.102 

0.221 

1 

o 

• 

l.. 5 

0.954 

5.242 

- 0.954 

0.191 

1 

o 

c_n 

1.5 

0.757 

5.003 

- 0.753 

0.151 

- 0.3 

1.5 

0.670 

4.570 

- 0.547 

0.110 

- 0.1 

1.5 

0.598 

4.034 

- 0.363 

0.0727 

0.1 

1.5 

0.547 

3.481 

- 0.209 

0.0420 

0.3 

1.5 

0.518 

2.956 

- 0.0852 

0.0171 

0.5 

1.5 

0.510 

2.465 

0.0171 

- 0.00343 

0.7 

1.5 

0.525 

1.981 

0.108 

- 0.0216 

0.9 

1.5 

0.572 

1 . 383 

0.212 

- 0.0426 

1.0 

1.5 


- 0 

-+■ CO 

->• - OO 

0.1 

1 . 1 

0.487 

3.377 

- 0.730 

0.146 

0.3 

1 . 3 

0.425 

2.978 

- 0.229 

0.0459 

0.7 

1 .7 

0.619 

1.961 

0.211 

- 0.0424 

0.9 

1 .9 

0.731 

1 .401 

0.41 2 

- 0.0826 

1.0 

2.0 


- 0 

-> 00 

- CO 











1 x 

v ( x , 0 ) = J S f(x)dx , 
a 2 


v(x,0) 


v 2 (x,+0) , 
v-| (x,+0) , 


f(x) 


f 2 (x) , 
f](x) , 


( a 2 1 x 1 b i ) > 

(a 2 < x < a) , 
(a < x < b ] ) , 
(a 2 < x < a) , 
( a £ x £ b-| ) . 


(25) 


The results given in this paper show the effect of the 
inclusion-crack geometry and the material constants on the be- 
havior of the stresses around the singular points. In addition 
to their application to fracture through conventional theories 
whenever valid, they may be used in connection with a simple 
criterion such as "a maximum cleavage strength at a character- 
istic distance" in studying fracture initiation from singular 
points where the power is not -0.5. It should also be noted 
that the problem of radial cracks which are not collinear may 
be solved without too much difficulty by using the technique 
described in this paper. 
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